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We report on the results of directional point-contact Andreev-reflection (PCAR) measurements 
in Ba(Fei_xCox)2As2 single crystals and epitaxial c-axis oriented films with x — 0.08 as well as 
in Ca(Fei_xCox)2As2 single crystals with x = 0.06. The PCAR spectra are analyzed within the 
two-band 3D version of the Blonder- Tinkham-Klapwijk model for Andreev reflection we recently 
developed, and that makes use of an analytical expression for the Fermi surface that mimics the 
one calculated within the density-functional theory (DFT). The spectra in Ca(Feo.94Coo.o6)2As2 
unambiguously demonstrate the presence of nodes or zeros in the small gap. In Ba(Feo.92Coo.o8)2As2, 
the a6-plane spectra in single crystals can be fitted by assuming two nodeless gaps, but this model 
fails to fit the c-axis ones in epitaxial films. All these results are discussed in comparison with recent 
theoretical predictions about the occurrence of accidental 3D nodes and the presence of "hot spots" 
in the gaps of 122 compounds. 



INTRODUCTION 

Four years after the discovery of superconductivity in 
Fe-based compounds [1] , a surprising variety of behaviors 
has been unveiled not only between compounds of differ- 
ent classes (the so-called 1111, 11 or 122 classes, just to 
mention the most studied) but also among compounds 
of the same class that share the same crystal structure. 
One of the properties that shows unexpected variations 
in a given class of Fe-based pnictides is the symmetry of 
the order parameter(s) (OPs) so that it is not possible 
to make a general statement about the gap symmetry 
of these compounds. The multiband electronic structure 
and the Fermi surface (FS) nesting that is responsible 
for the antiferromagnetic instability in the parent com- 
pounds strongly suggest a superconducting pairing me- 
diated by spin fluctuations [2|,l3l which is also supported 
by many experimental facts [4|. Various possible gap 
symmetries are allowed, at least theoretically, within the 
picture of spin-fluctuation mediated pairing. In most of 
optimally-doped compounds the so-called 5± symmetry 
with isotropic OPs of opposite sign on the holelike and 
electronlike FS sheets is realized [2], but the emergence 
of nodes or zeros on some FSs is theoretically predicted 
in various situations [5-7]. 

Here we report on recent PCAR spectroscopy results in 
single crystals of 8% Co-doped BaFe2As2 (Ba-122) with 
bulk T^^^ = 24.5 K and in single crystals of 6 % Co-doped 
CaFe2As2 (Ca-122) with T^^^ = 20.0 K. In the first case, 
the PCAR spectra (taken with the current injected along 
the ah planes) are compatible with two nodeless gaps of 
different amplitude, the smaller being probably associ- 
ated to one of the electrolike FS sheets, in agreement with 



angle-resolved photoemission spectroscopy ( ARPES) j8| . 
In the crystals of Co-doped Ca-122, the PCAR spectra 
systematically show a zero-bias peak or maximum, which 
strongly suggests that the small gap (residing here on the 
holelike FS sheet centered at the F point of the Brillouin 
zone) presents lines of nodes or zeros on the FS, while 
the large gap remains presumably isotropic. This finding 
agrees with the predictions about the emergence of 3D 
nodes in the OP of 122 compounds when the holelike FS 
evolves towards a topological transition from a warped 
cylinder to separate pockets 0, [sl • Finally, we present 
new results of PCAR experiments in epitaxial thin films 
of 8% Co-doped Ba-122 (with TJ^^^ = 24.35 K). In this 
case, the current is injected along the c axis and the 3D- 
BTK model with two nodeless gaps fails to fit the PCAR 
spectra at low energy. This might be a hallmark of the 
presence of "hot spots" (where the gap is strongly sup- 
pressed) or 3D nodal lines on some of the FSs, which 
has been proposed [5] to explain Raman results in this 
compound [10|]. 



EXPERIMENTAL DETAILS 
The samples 

The Ba(Fei_xCox)2As2 single crystals (with x = 0.08) 
were prepared by the self- flux method [11] under a pres- 
sure of 280 MPa at the National High Magnetic Field 
Laboratory in Tallahassee. The crystal size is ~ 1x1x0.1 
mm^, and the c axis is perpendicular to the larger sur- 
face. The resistive transition sets in at T^^ = 24.5 K 
with ATc (10%-90%) = 1 K. 
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The Ba(Fei_xCox)2As2 films (with x = 0.08) were 
deposited at the Leibniz Institute for Solid State and 
Materials Research (IFW) in Dresden, Germany. Two 
kinds of substrates were used: single-crystalline CaF2 or 
MgO (in this case, with a Fe buffer layer on top of it). 
The films are epitaxial with the c axis normal to the 
surface. The superconducting transition of the films on 
CaF2 has a midpoint at T^'"^ = 24.35 K, while its width 
is ATc = T^^o% - T^io% = 1.7 K. For the films on MgO 
the same quantities are TJ^*^ = 23.80 K and ATc = 1-5 
K. 

The Ca(Fei_a;Coa^)2As2 single crystals were grown at 
ETH Zurich from Sn flux, as described in Ref. p^], and 
were plate- like, with the c axis perpendicular to the plate. 
The superconducting transition as measured from DC 
susceptibility sets in at T^^ = 20 K and an effective 
T^ff = 17 K can be determined by extrapolating the 
linear part of the curve. The transition is rather broad, 
but this is common to all the state-of-the-art crystals of 
the same compound grown at present [l^. The broad 
transition is however not detrimental to point-contact 
measurements: Simply, the local critical temperature of 
point contacts can vary from point to point within the 
transition width. 



Point-contact setup and procedure 

The point contacts were made by using the "soft" tech- 
nique described elsewhere [l^, in which the metallic tip 
used as counterelectrode in the more conventional needle- 
anvil configuration is replaced by a small drop of Ag paste 
{(j) :^ 50 /im) put on the sample surface. In single crystals 
the contacts were put on a fresh side surface (exposed by 
breaking the crystals) so that the current was mainly in- 
jected along the ab planes. In films, the point contacts 
were put on the top surface. Owing to the film orien- 
tation, this means that the probe current was mainly 
injected along the c axis. 

The conductance curves d//dV vs. V of each point 
contact were obtained by numerically differentiating the 
measured I — V characteristics. The spectra were then 
normalized (i.e. divided by the normal-state conductance 
of the same contact) to allow a comparison to the rele- 
vant theoretical models. In single crystals (where the 
high-bias tails of the spectra at different temperatures 
are superimposed to one another) we divided by the con- 
ductance curve recorded just above the critical temper- 
ature of the contact, (defined as the temperature at 
which the Andreev reflection features disappear). In the 
case of films, this procedure is impossible because of an 
abrupt vertical shift of the spectra when the temperature 
approaches T^. This shift is related to the transition 
to the normal state of the portion of the film between 
the point contact and the voltage electrode [l5| and has 
nothing to do with the contact itself (a detailed model 



of this phenomenon will be given in a forthcoming pa- 
per). Therefore, the low-temperature experimental spec- 
tra were divided by a guess normal state obtained by 
fitting their high-energy tails with a fourth-order poly- 
nomial function, as we did in Ref.[l6[. 

CALCULATION OF THE FERMI SURFACE 

The FS of the materials under study 
(Ba(Fei_xCox)2As2 with x = 0.08 and 
Ca(Fei_xCox)2As2 with x = 0.06) was calculated 
within the Density-Functional Theory (DFT) by using 
the Elk FP-LAPW Code (http : //elk.sourceforge.net/) 
and the GGA approach for the exchange correlation 
potential. To account for the partial substitution of Fe 
with Co we used a virtual-crystal approximation that 
indeed works particularly well for the Fe-Co substitution. 

The lattice constants used for the calculations in the 
case of Ba(Fei_xCox)2As2 are a = b = 3.9625 A and 
c = 13.0168 A as in Ref. [17]. The height of the As atom 
above the Fe layer is Has = 1.1975 A. The resulting FS 
is shown in Fig. [TJa). It features two holelike FS sheets 
at the center of the Brillouin zone (BZ). Both have the 
shape of warped cylinders whose cross section is maxi- 
mum at the top and bottom edges of the BZ, but the 
outer one shows a more marked warping. The situation 
is somewhat similar to that found in BaFe2As2 with par- 
tial substitution of As with P, as described by Suzuki et 
al. The two electronlike FS sheets at the corners of 
the BZ are also warped cylinders with the characteristic 
elliptical cross section whose semimajor axis varies along 

In the case of Ca(Fei_xCox)2As2, there is no direct 
experimental information on the low-temperature lattice 
constants. Owing to the small dependence of the room- 
temperature lattice parameters on the doping content 
[18i] , we assumed the low-temperature lattice constants of 
the parent compound CaFe2 A2 in the tetragonal phase to 
be a good first approximation to the real ones at the dop- 
ing content of our interest {x = 0.06). We then started 
from the lattice constants of the orthorhombic phase of 
CaFe2A2 at pressure P = calculated as in Ref. [l9| , 
and then we made the structure tetragonal by averaging 
a and b. The result is a = b = 3.925 A and c = 11.356 A. 
These values are in good agreement with the experimen- 
tal ones measured in the tetragonal phase of CaFe2A2 
at 300 K (and pressure P = 0.8 - 1 GPa [20|). Starting 
from the calculated equilibrium phase and always consid- 
ering the antiferromagnetic phase, an optimized param- 
eter Has = 1.309 A was obtained. The charge density 
was thus integrated over 8 x 8 x 4 /c points in the Bril- 
louin zone and the band structure as well as the FSs were 
calculated in the non-magnetic body-centered tetragonal 
phase. The resulting FS is shown in Fig. [U^c). It is clear 
that at this doping content the holelike FS sheets are un- 
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dergoing a topological transition. While at lower doping 
they have the shape of strongly warped cylinders (similar 
to those shown for Ba-122 in Fig. [U^a)), at x = 0.06 they 
split into separate cup-shaped pockets centered around 
the Z points. 



RESULTS AND ANALYSIS 

Ba(Fei_ccCocc)2As2 single crystals 

Figure [2] presents two examples of experimental nor- 
malized conductance curves (symbols) of a6-plane con- 
tacts in Ba(Fei_xCox)2As2 single crystals. Both the 
curves in (a) and (b) show a zero-bias dip, two sym- 
metric maxima related to a small gap, and two other 
kinks related to a second, larger gap. Moreover, the 
spectra feature additional, higher-energy structures (in- 
dicated by arrows) that are likely to be due to the strong 
electron-boson coupling [2l| . We will not deal with these 
structures here; they cannot be reproduced by any BTK 
model unless an energy-dependent gap (as obtained from 
the solution of Eliashberg equations) is used instead of 
the constant BCS-like gap used here [2l|, [22|. To fit the 
experimental curves in Fig. [2] we used the 3D version 
of the generalized BTK model j22[, that accounts for 
the real shape of the FS and is not based on the sim- 
plifying assumption that the latter is spherical or cylin- 
drical, as instead the ID and 2D BTK models do. For 
ease of calculation, we used the model FS shown in Fig. 




FIG. 1: (a) The Fermi surface (FS) of Ba(Fei-xCox)2As2 
with X = 0.08. (b) The model FS used in the 3D-BTK model 
to fit the experimental conductance curves. Solid surfaces 
represent the FS sheets, while gridded surfaces represent the 
amplitude of the corresponding energy gap. (c) The FS of 
Ca(Fei_xCox)2As2 with x = 0.06. (d) The model FS (in the 
upper half of the first Brillouin zone) used in the 3D-BTK 
model to fit the experimental conductance curves. Gridded 
surfaces represent the amplitudes of the gaps, a large isotropic 
one on the electron FS and a small anisotropic one on the hole 
FS. 



[TJb), which mimics the real one from DFT calculations 
reported in Fig. [IJa). The former consists of two sep- 
arate hyperboloids of revolution, meant to simulate the 
main holelike and electronlike sheets, whose radii at the 
center and at the top (and bottom) of the BZ are in the 
same proportions as in the real FS, although in Fig. [Hb) 
the distance between them has been enhanced for clar- 
ity. The normalized conductance expected in the case 
of current injection along the ah plane can be calculated 
by using eq. 9 of Ref. [22]. Note that, due to the par- 
ticular mathematical form of the equation, the relative 
position of the two FS sheets is irrelevant. Of course, 
the choice of the symmetry of the order parameters that 
reside on the two hyperboloids is determined by the ex- 
periment. In this case, the absence of zero-bias maxima 
(in all the spectra we measured) suggests that both the 
gaps are likely to be nodeless, in agreement with ARPES 
measurements [8]. We therefore used in our model two 
5-wave gaps whose amplitudes are pictorially indicated 
as gridded surfaces in Fig. [IJb). We will call Ai and A2 
the large and the small gap, respectively. The 3D BTK 
model contains 3 adjustable parameters for each band, 
i.e. the gap amplitude A, the broadening parameter F 
and the barrier parameter Z . The weight of each band 
in the conductance is not adjustable (as instead happens 
in the 2D model) being automatically determined by the 
geometry of the FS and by the barrier parameter Z. This 
considerably reduces the uncertainty on the best-fitting 
parameters. 

Fig. [2] shows two examples of fit (lines) to the ex- 
perimental spectra. The values of the parameters are 
indicated in the caption; the average values are Ai = 
11.75±0.25 meV and A2 = 3.8±0.8. In the 2D-BTK fit of 
the same curves the gaps turned out to be Ai = 10.7±0.2 
meV and A2 = 4.4 ± 0.6 meV. With respect to the 2D 
BTK fit, the 3D one gives smaller values of the small gap 
and larger values of the large gap, and the theoretical 
curves are narrower. This is due to the fact that in the 3D 
model the weight is fixed (in these two cases, the weight 
of band 1 is about 0.2). Indeed, the 2D fit can be forced 
to follow the experimental curve at energies higher than 
13 meV (where the 3D fit fails) if the weight of the bands 
is kept around 0.5, but this clearly would not reflect the 
real shape of the FS sheets. The inability of the model 
to reproduce the higher-energy structures (in particular 
the kinks at about 20 meV) simply confirms that these 
structures are related to effects that are not accounted 
for by the model (and indeed can be explained as being 
due to the strong electron-boson coupling, as shown in 
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Ca(Fei_a:Cox)2As2 single crystals 

Figure [3] shows two examples of typical conductance 
curves of a6-plane contacts in Ca(Fei_xCox)2As2 single 
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crystals. Contrary to what happens in Co-doped Ba-122, 
here 100% of the conductance curves presents zero-bias 
maxima or peaks. As shown elsewhere [22|, this is a clear 
sign that one of the gaps is strongly anisotropic in the 
(kxky) plane. According to theoretical calculations 0, 
nodal lines can appear in the order parameter of the outer 
holelike FS sheet when the size of the latter is largely 
increased in the vicinity of the top and bottom faces of 
the BZ. In Ref. this effect is produced by the reduction 
in the pnictogen height Has induced by the substitution 
of As with P. Our DFT calculations show that the effect 
of Co substitution in CaFe2As2 is very similar, and that 
the FS of Fig. [TJc) is actually the extreme consequence of 
a doping-induced increase in the warping of the holelike 
sheet. Thus, our PCAR measurements confirm that, even 
within a general 5± picture of spin-fluctuation mediated 
superconductivity, nodal lines can appear in the holelike 
FS when the latter is strongly deformed - in particular, 
if it undergoes a topological transition and splits into 
separate closed pockets. 



We thus modeled the FS with one hyperboloid (for 
the electronlike sheets) and one spheroid (for the hole- 
like pocket) as in Fig. [TJd) and assumed an isotropic 
large gap Ai on the former and an anisotropic small gap 
A2 on the latter. This gap should represent an evolution 
of that depicted in Refs. and [23| when the relevant 
FS sheet becomes a closed surface. One possibility is 
to use a (i-wave gap A2(6>,0) = A2 cos(26>) sin(^) as in 
Ref. [sl, to reproduce the sign change of the OP. How- 
ever, the symmetry of the OP of Refs. and [23| in 
the {kxky) plane is different from the (i-wave one and 
features two additional axes of equation ky = i/c^; like 
the s + Qx'^-y'^ oiie [24]. As a consequence, the prob- 
ability of constructive interference between electronlike 
and holelike quasiparticles (that gives rise to the zero- 
bias peak) is strongly reduced. However, the zero-bias 
maximum could also arise from the existence of angular 
regions where the gap has very small amplitude, rather 



than from its change of sign, as shown in Ref. [22|. To 



account for this possibility, here we modeled the gap on 
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FIG. 2: Two examples of normalized conductance curves of 
a5-plane point contacts in Ba(Fei-xCox)2As2 single crystals 
with X — 0.08 (symbols) and of the relevant fit (lines) with 
the 3D BTK model using the FS shown in Fig. ^h). The 
fitting parameters of the curve in (a) are: Ai = 11.5 meV, 
Fi 1.85 meV, Zi 0.03 and A2 3.0 meV, F2 3.6 meV, 
Z2 — 0.310. The weights of the two bands in the conductance 
are wi = 0.23 and W2 = 0.77. The parameters for the fit of 
the curve in (b) are instead Ai = 12.0 meV, Fi = 1.75 meV, 
Zi 0.08 and A2 4.6 meV, F2 3.00 meV, Z2 0.245. 
The weights in this case are wi = 0.21 and W2 = 0.79. 



FIG. 3: Two examples of normalized conductance curves of 
a6-plane point contacts in Ca(Fei-xCox)2As2 single crystals 
with X = 0.06 (symbols) and of the relevant fit (lines) with 
the 3D BTK model using the FS shown in Fig. ^d). The 
fitting parameters of the curve in (a) are: Ai = 5.7 meV, 
Fi = 6.15 meV, Zi = 0.145 and A2 = 1.4 meV, F2 = 1.25 
meV, Z2 — 0.05. The same parameters for the fit of the curve 
in (b) are instead Ai = 5.6 meV, Fi = 6.5 meV, Zi = 0.23 
and A2 = 1.45 meV, F2 = 1.7 meV, Z2 = 0.05. In both (a) 
and (b) the angle between the normal to the interface and the 
a axis is a = 7r/8. 
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the holelike pockets with a fully anisotropic 5-wave gap 
of equation A2(^,^) = A2 cos^(2^) sin(^) (gridded sur- 
face in Fig. md)). The fit of the experimental curves is 
shown in Fig. [3] (solid lines). The parameters of these 
two particular fits are listed in the caption. From differ- 
ent fits we get the following average values for the gaps: 
A2 = 1.4 ± 0.1 meV and Ai = 5.5 ± 0.3 meV. If one 
chooses a d-wave symmetry for the small gap, a fit of 
comparable quality is obtained but the values of the gaps 
are Ai = 1.6 ± 0.1 meV and A2 = 5.3 ± 0.2 meV 



Ba(Fei_a,Cox)2As2 films 

Two examples of PCAR spectra in Ba(Fei_xCox)2As2 
films are shown in Fig. HI The spectrum in (a) refers to 
a film on MgO substrate with Fe buffer layer, and shows 
double-gap features similar to those of Fig. [2jb) but 
even clearer. The spectrum of panel (b) was measured in 
a film on CaF2 and presents much smoother structures. 
In particular, the large gap manifests itself in very broad 
shoulders at about ±6 meV. As stated in a previous sec- 
tion, the normalization in this case gives rise to some 
ambiguity since the normal state at is not related to 
the contact alone, but includes a contribution from the 
portion of the film between the point contact and the 
voltage electrode. The normalization by a guess normal 
state can be used, but one should keep in mind that the 
choice of the presumed normal state can affect the shape 
of the curve and, to some extent, the fitting parameters 
- actually, the broadening F and the barrier parameter 
Z are the most affected by that choice while the values 
of the gaps are rather robust. The dashed lines in Fig. |4] 
represent the results of the 3D BTK fit using the same 
model FS as in Fig. [TJa). The 3D model is clearly unable 
to fit the curve very well because of a Z-enhancing effect 
due to the shape of the FS [22| . The parameters of the fit 
shown in Fig. H^a) are Ai = 10.8 meV, Fi = 3.85 meV, 
Zi = 0.1, A2 = 3.0 meV, r2 = 1.4 meV, Z2 = 0. The 
weights of the bands, which are not adjustable within this 
model, are wi = 0.39 and W2 = 0.61. Similarly, in Fig. 
m^b) the 3D BTK fit (dashed curve) was obtained by us- 
ing Ai = 7.3 meV, Fi = 2.65 meV, Zi = 0.02, A2 = 3.3 
meV, r2 = 1.2 meV, Z2 = 0. Note that in both cases a 
very deep zero-bias minimum is obtained despite the fact 
that the intrinsic barrier strength is put equal to zero in 
band 2 and is very small in band 1. The failure of the 3D 
BTK model may be due to different factors. One possibil- 
ity is that, for some reasons, the probe current is injected 
mainly along the ab plane despite the orientation of the 
film. This may happen, for example, when the surface is 
rough, a minority percentage of grains are not oriented 
or the point contact is made by using a sharp tip that 
may pierce the surface. In our case, all these possibilities 
are rather unlikely. The other possibility is that the use 
of two isotropic gaps is not the best choice to reproduce 



the real distribution of the gaps over the FS. Indeed, the 
presence of "hot spots" where the gap is significantly sup- 
pressed has been evidenced by Raman spectroscopy [10| 
and then justified theoretically [5]. Since the hot spots 
seem to reside on the electronlike FS sheets, one might 
try to improve the model by using a small anisotropic 
gap with zeros and a large isotropic gap (as we did in 
the case of Ca(Fei_xCox)2As2). This would certainly 
improve the fit because there would be quasiparticle ex- 
citations even at very low energy and the Z-enhancing 
effect might be compensated by the zero-bias maximum 
due to these excitations [l^. We will explore this possi- 
bility elsewhere. For the time being, we show in Fig. 2] 
what would be the result of a 2D-BTK fit (solid lines). 
The fit is rather good in the energy region of the gaps 
and catches the main features of the curves, but just be- 
cause the weights of the bands were taken to be different 
from those determined by the shape of the FS in the 3D 
fit. In (a) we got Ai = 9.5 meV and A2 = 5.5 meV by 
using wi = 0.44 and W2 = 0.56. Of course, changing the 
weight affects all the fitting parameters: the resulting gap 




o 

B 
o 

"D 

o 
o 

"D 
CD 
_N 

IB 

E 
o 




-40 -30 -20 -10 10 20 30 40 










: 7.1 meV 




4.43 meV 


^ = 


: 7.3 meV 




3.3 meV 



£ 1.17 



1.08 



0.99 

-40 -30 -20 -10 10 20 30 40 

Voltage (mV) 

FIG. 4: Two examples of normalized conductance curves of c- 
axis point contacts in Ba(Fei-xCox)2As2 films with x = 0.08 
(symbols) but on different substrates: MgO + Fe buffer layer 
in (a) and CaF2 with no buffer layer in (b). The dashed lines 
represent the 3D BTK fit, while solid lines represent the 2D 
BTK fit. The corresponding gap amplitudes are indicated 
in the labels. The weights of the bands in the 3D fit are 
wi = 0.39 and W2 = 0.61. The 2D fit was performed by fixing 
instead wi — 0.44 and W2 — 0.56. 
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values (averaged over different fits) are Ai = 9.5 ±0.2 
meV and A2 = 5.4 ± 0.2 meV. The uncertainty is small 
because the peaks are rather sharp. In (b) we obtained 
Ai = 7.1 meV and A2 = 4.43 meV, with the weights 
wi = W2 = 0.5. Averaging over different fits obtained 
by changing the weight, the gap values turn out to be 
Ai = 7.4 ± 0.3 meV and A2 = 4.2 ± 0.25. Note that 
the 2D-BTK fit of similar curves (with current injection 
along the c axis) in single crystals of the same material 
gave Ai = 9.2 ±1.0 meV and A2 = 4.1 ± 0.4 meV [21|]. 
There is thus a partial superposition of the gap values, 
although not a perfect agreement. This is not particu- 
larly surprising since: i) there may be some effect of the 
substrate, especially in the case of the films on Fe buffer 
layer; ii) the measurements in films highlight a certain 
degree of inhomogeneity in the superconducting proper- 
ties, although the film itself is of very high quality; iii) 
the statistics of the measurements in films needs to be 
extended to draw definite conclusions. 



fit of the c-axis spectra measured in thin films with the 
same model is unsuccessful around zero bias, where the 
so-called Z-enhancing effect due to the shape of the FS 
produces a strong depression of the conductance which is 
not observed experimentally. This might be the sign of 
the presence of gap minima ^ whose existence has been 
inferred from Raman spectroscopy results [lo|. 
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CONCLUSIONS 

We have presented the results of PCAR experiments 
in Ba(Fei_xCox)2As2 single crystals and films with x = 
0.08 and in Ca(Fei_xCox)2As2 single crystals with x = 
0.06. The results of these measurements have been an- 
alyzed by using a novel 3D version of the BTK model 
that accounts for the shape of the FS - actually using an 
analytical model for the real FS calculated within DFT. 
The results indicate that the two systems, although be- 
longing to the same 122 family of Fe-based compounds, 
show major differences in the shape of the spectra, and 
this is a symptom of a difference in the symmetry of the 
order parameters. In Ca(Fei_xCox)2 As2 a zero-bias peak 
or maximum is always observed and can be explained as 
being due to a strong anisotropy of the smaller gap. DFT 
calculations show that, at 8% Co content, the holelike 
FS sheet undergoes a topological transition and splits 
into separated closed pockets. This can be interpreted 
as an extreme consequence of the increasing- with-doping 
warping of this FS sheet, which is predicted to be ac- 
companied by the emergence of 3D nodes in the relevant 
gap. Indeed, a fit of the spectra is possible by assum- 
ing either a d-wave [9] or an anisotropic symmetry of 
the small gap, to mimic the possible evolution of these 
3D nodes when the topological transition occurs. In the 
case of Ba(Fei_xCox)2As2 at x = 0.08, the calculated 
FS shows markedly warped holelike sheets but no topo- 
logical transitions. The spectra taken in single crystals 
with the current injected along the ab plane can be fitted 
with two isotropic gaps; the agreement between model 
and data is limited to the central region of the spectra 
because of additional structures around 20 meV that are 
due to the strong electron-boson coupling 
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